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The well known interpretational difficulties with nonlinear Schrodinger and von Neumann equa- 
tions can be reduced to the problem of computing multiple-time correlation functions in the absence 
of Heisenberg picture. Having no Heisenberg picture one often resorts to Zeno-type reasoning which 
explicitly involves the projection postulate as a means of computing conditional and joint probabili- 
ties. Although the method works well in linear quantum mechanics, it completely fails for nonlinear 
evolutions. We propose an alternative way of performing the same task in linear quantum mechanics 
and show that the method smoothly extends to the nonlinear domain. The trick is to use appropriate 
I time-dependent Hamiltonians which involve "switching-off functions". We apply the technique to 

■ the EPR problem in nonlinear quantum mechanics and show that paradoxes of Gisin and Polchinski 

(N . 

disappear. 
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I. INTRODUCTION 



Ten years ago Polchinski published a paper jl] where he showed how to avoid unphysical influences between 
separated systems in nonlinear quantum mechanics. The paper was a response to criticisms of Kibble- Weinberg 
nonhnear quantum mechanics |^,^ raised by Gisin [Q and one of us |^] . Polchinski proposed to eliminate the problem 
by eliminating its source, the projection postulate, and resorting to the Many Worlds Interpretation. However, after 
having eliminated one difficulty he showed that there appears another one, a kind of weird communication between 
^^O . different Everett worlds. 

' The paper had several implications. First of all, it convinced a large group of physicists working on foundations of 
quantum mechanics, including Weinberg himself (sec his remarks in [^), that linearity is an ultimate element of all 
the possible "final theories" . 

I On the other hand, the Polchinski multiparticle extension showed that at least part of difficultes trivally disappears if 
one starts with nonlinear evolutions formulated in terms of von Neumann-type equations. Looking at the von Neumann 
dynamics as a classical Hamiltonian system on a Poisson manifold of states represented by reduced density matrices, 
one arrives at a surprisingly elegant formalism [|7|-|l2| which goes much beyond the original Weinberg proposal. Theories 
^2 ' so constructed led to new families of integrable von Neumann equations which turned out to be very interesting in 
qh, themselves at least from the point of view of integrable systems [|l3|-0 . 

In spite of the original problems with the Weinberg formalism notable developments were reported by groups working 
on "geometric quantum mechanics" JTsj-pOt. Quantum mechanics, when looked at from a geometric standpoint, turns 
out to be a classical theory "living" on a very symmetric phase space (a Kahler manifold with maximal symmetries). 
5— i ' It becomes quite natural, then, to think of standard quantum mechanics as an analogue of a "vacuum" spacetime from 
general relativity. Since the most interesting spacetimes are those which are less symmetric, maybe the same will be 
true for quantum mechanics? An unquestionable beauty of the geometric perspective suggests that these possibilities 
should be further investigated. 

Yet another road to nonlinear Schrodinger equations was discovered by Goldin and Doebner pl]-|2^. Studying 
representations of groups of diffeomorphisms they noticed that there exists a class of representations parametrized 
by the real number D which has an interpretation of a diffusion coefflcient |^,^. All the other representations had 
a clear physical interpretation. They were known to unify the description of an extraordinary variety of quantum 
systems and generated some physical predictions difficult to reach by other means (e.g. the anyon statistics in 
two-dimensional space). However, this particular class resisted a quantum mechanical interpretation until it was 
realized that it corresponds to a family of nonlinear dissipative Schrodinger equations earlier derived in the context of 
quantum chemistry by Schuch and his collaborators [ pTj . Doebner and Goldin subsequently introduced the group of 
nonlinear gauge transformations necessary to understand the resulting theory. The resulting formalism was shown to 
contain as special cases nonlinear modifications of quantum mechanics proposed independently by other researchers 
p8|-p4| It became clear that linearity is a gauge-dependent property and a simple rule "linear good, nonlinear bad" 
is oversimplified. 

Quite recently it turned out that dynamics of certain Z?-branes leads to nonlinear Schrodinger equations od the 
Doebner-Goldin type [pa,p6| and the von Neumann-type nonlinearities introduced in have links with nonextensive 



1 



statistics IfQ. 

In spite of the optimism of the groups we have mentioned not everyone was fully satisfied by the Many Worlds 
solution. Gisin openly disqualified nonlinear Schrodinger equations as "irrelevant" |4^. Goldin |^l| did not agree 
with Gisin but formulated a series of doubts with respect to the Polchinski solution pointing out, in particular, that 
multiparticle extensions of physically equivalent theories but performed in different nonlinear gauges lead to physically 
inequivalent results. Also Mielnik in a recent paper argued on general grounds that the difficulties remain. 

As a result, a person who wants to understand the current state of knowledge about nonlinear generalizations of 
quantum mechanics will find a complete spectrum of mutually contradicting statements. 

We think that three main questions are causing the entire embarassment. One is the unclear role of nonlinear aver- 
ages and eigenvalues if nonlinear observables are involved. The other is the role of density matrices ( "mixed" states) 
in nonlinear formalism. Finally, an open question is how to compute joint probabilities in EPR-type experiments. 

The main topic of the present paper is the third issue. Two applications of nonlinear averages are discussed in the 
Appendix. We do not address the problem of density matrices. We believe that a common error one finds in the 
literature on nonlinear Schrodinger equations is to treat solutions of von Neumann-type equations as states which are 
mixed in the standard sense of mixed states on a phase space ji^ . To the best of our knowledge no approach based 
on Schrodinger-type equations was able to successfully address the question of composite systems. This should be 
contrasted with the formalism based on von Neumann equations which, up to a few interpretational subtleties which 
are a subject of the present work, simply does not lead to this difficulty. 

In what follows we formulate a new operational framework for interpretation of experiments in both linear and 
nonlinear quantum systems. We do not use the Many Worlds Interpretation and hence do not have the "Everett 
phone" . On the other hand we make a very restrictive use of the projection postulate. We allow to use projections 
only in two cases: (1) to determine initial conditions and (2) at the moments particles are detected and destroyed 
(final conditions). 

The question that arises, then, is how to describe situations where there are, say, two correlated particles and one 
of them is destroyed by a detector at time t — ti whereas the other one is allowed to evolve until it gets destroyed at 
t = t2 > ti. This is what we call a two-time measurement. A typical solution, known also as a quantum Zeno effect, 
is to regard the second particle as evolving for ti < t < t2 with a new initial condition produced by an appropriate 
projection-at-a-distance at t — ti. 

We exclude this approach since in application to nonlinear evolutions of entangled states it leads to serious compli- 
cations. Instead, we reformulate the standard strategy with the help of "switching-off functions" and show that (1) 
in the linear case the result is identical to the standard Zeno-type calculation and (2) in nonlinear case is free from 
unphysical influences between correlated particles and, hence, ultimately eliminates the "EPR phone" . 

We will begin with a few remarks on filters and histories. The notion of a filter often occurs in the context 
of conditional probabilities and the projection postulate. However, there exists also a desription where filters are 
described in terms of unitary matrices and do not lead to irreversible projections. This point is essential for our 
argument. In the next two sections we discuss the role of the projection postulate for two-time measurements and 
point out the role of "switching-off functions" . In Sec. V we discuss problems with two-time measurements in presence 
of nonlinearities. Sec. VI discusses Polchinski's multi-particle extension and its modification in terms of switching-off 
functions. In Sec. VII we explicitly solve an example and compare the standard projection-based calculation with 
the new proposal. Sec. VIII discusses implications of our approach for the question of teleportation. Arguments are 
given that one should not identify pre- and post-selected ensembles if nonlinear evolutions are involved. The point is 
essential since a true teleportation corresponds to a pre-selected ensemble whereas ordinary correlation experiments 
involve post-selections. 

We end with the Appendix on nonlinear averages, the use of them in information theory and nonextensive statistics, 
and their links with the subject of the paper. 

II. HISTORIES AND FILTERS 

A typical quantum history is shown in Fig. 1. A particle produced by a source undergoes a series of scattering 
events, represented by unitary evolutions Ui and U2- Before, between and after the scattering the evolution is free, 
say. The particle may, for simplicity, be regarded as a two level system and the Us as some sorts of beam splitters 
(mirrors, analyzers of polarizaton, Stern-Gerlach magnets,...). There are several outgoing channels, and each of the 
channels corresponds to a different result of the experiment. 

The probability of detection in, say, channel 1 can be calculated in a couple of different ways. One way of doing 
this is to compute the overall evolution operator V{t), including the whole evolution from the source till detection, 
and calculate the probability 
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Pi = Tr {E2Vit)p{0)V{ty) 
= Tr {E2p{t)) 



(1) 



where E2 is an appropriate projector operator. 

Another way is to spht the evolution into pieces exphcitly involving the scattering matrices Uk and a free evolution 
between them. The Us play here a role of analyzers of the properties represented by projection operators Ek and 
'^Ek, Ek + ^Ek — I ■ Denote by U{t — t') the free evolution between times t and t' . Since it satisfies the group 
composition property 

U{t-t')^U{t')^U{t). (2) 

we find 

Pi = Tr (U{t - t2)E2U{t2 - ti)EiU{ti)p{Q)U{ti)^EiU{t2 - ti)^E2U{t - ta)^) (3) 

= Tr(E2{t2)Ei{ti)p{Q)Ei{ti)E2{t2)) (4) 

where Ek{tk) — U {t^)'^ EkU [t]^) . Expression (||) is often called a history In quantum-optics literature one would 
rather speak of a two-time correlation function. In derivation of (^) one makes an explicit use of the unitarity of 
evolution whereas the transition from (^) to (Q) involves (||). On the other hand, the form (|l]) is more general and 
can survive even if the dynamics 1 1— > p(t) is nonlinear. 

It is an appropriate moment to point out that one should not treat the first splitting of the beam by Ui as a 
measurement of the property Ei . The measurement occurs when the particle vanishes in the detector and produces 
a click. Before that moment the evolution is unitary and reversible. After the measurement the measured particle is 
destroyed and is no longer under the jurisdiction of quantum mechanics. 

Our beam splitters Uk are examples of devices which are sometimes referred to as filters |4^] . As we have seen, the 
derivation of the history explicitly employed the fact that, for calculational purposes, filters can be represented 
by projectors and, hence, by non-unitary operators. However, it is extremely important to understand that the same 
physical process can be described by the exact formula which describes propagation of particles through filters in 
a unitary way. 

In such a formulation projectors are applied only twice — at the very beginning when one selects a concrete initial 
condition, and at the very end where the particle is detected and destroyed. 



III. ENTANGLED STATES AND TWO-TIME MEASUREMENTS 



Two-time measurements are at the heart of difficulties with nonlinear extensions of quantum mechanics. We begin 
with their orthodox description in linear quantum mechanics. 

Consider a two-electron system initially prepared in the singlet state 

|vI/(0))=.4sf|+)|-)-|-)|+)) (5) 

and whose Hamiltonian is 



The dynamics of the state is 



H = Hi (g) I2 + h (8) H2. (6) 

\^it)) = v,{t)^v2{t)\^m. (7) 



where Vk{t) = exp{—iHkt). It is known that |^'(0)) will be unchanged if Vi = V2. However, if Hi 7^ H2 the perfect 
initial anti-correlation of spins holding at t = will be lost at later times. 

Now let us imagine that two measurements of the z-components of spin are performed; one at time t = ti on particle 
#1, and the other ai t = t2 > ti on particle #2. We want to calculate an average of the observable az 'S) dz- 

Denoting by E± ® I and / ® E± the appropriate projectors we can proceed as follows. 



At t = ti project 



E±(8l\^{tr)) 
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• Evolve the resulting state by 

E±<»V2it2~h)\^{ti)) 

II E±(g>I\-^{h)) II 



(9) 



• Calculate at i = i2 the (conditional) probabilities 

{^{h)\E± V2{t2 - hyE±,V2it2 - h)\^{ti)) 

{^>{t,)\E±(g>I\^{t,)) 



(10) 



The denominator in ( |l0| ) is the probability of finding ±1 for the first particle at t = ti. Therefore the joint 
probability of finding ± for the first particle and ±' for the second one is the numerator of (p^), i.e. 

{^iO)\Viih)^E±Vi{h) ® V2it2)^E±,V2{t2)\^m. (11) 



The second step (^ was done in this way since we knew already the result of one measurement after (^). The result 
created an initial condition for the remaining dynamics of the other particle at i = ti (this is essentially the quantum 
Zeno effect). The projection postulate has therefore explicitly been used in the second step at i = ti. 

The formula ( pT] ) is again a correlation function. Let us note that we would get the same formula if we computed 
the average of the projector 

E± (g) E±' (12) 

in the state 



Vi{h)<»v2{t2)\-^m. (13) 

An identical result would have been obtained if one considered the time dependent Hamiltonian 

Ht,^tAt) = S{t-ti)Hi®l2 + e{t-t2)h®H2. (14) 
were the 6{x)s is the step function equal 1 for x <0 and otherwise, and computed the average of (|l^) in the state 

\^t,M{-t)) = g-'^'^^'^/o ^("-*^)'^"-^^^^^^/o «("-*^''^"|^(0)). (15) 

The times ti and t2 play a role of parameters describing the system. They naturally incorporate into the formalism 
the presence of detectors and the fact that detections are registered at certain times. 

It is easy to understand why the energy of the particles is not conserved during the dynamics: At the moments 
the particles "get frozen" by the detection process their energies are transferred to the detectors. Only the composite 
system "particles plus detectors" is closed. 

It is also obvious that one should not apply "switching-off functions" 6 for acts of filtering that do not end up in a 
detection and destruction, such as those denoted by Ui and Ei at Fig. 1. 



IV. TWO-TIME MEASUREMENTS ON BEAMS OF PAIRS 



In real experiments one deals with beams of pairs. A beam consisting of N pairs is a tensor product of N copies 
of single-pair states. Real experiments do not produce all the pairs at the same initial i = 0. An experiment may 
run for weeks, therefore it is much more realistic to impose initial conditions on each of the pairs separetely and at 
different times: tj, . . . , t^. 

A run consisting of N pairs of detections at times {t\^t\)^ ... (t]^,t^), can be described by the state vector 

If experimental conditions do not change during the run one may parametrize (p^ ) by times of flight At^ = ^fc ~ ^o- 
Observables, such as the one from the previous section are now represented in a frequency form i.e. 
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FNidz (g) cr^) = — {^{a^ cr^) (g) (J® I)® ... (g) (/«)/) + ... + (/«)/)«)... (g) (/ (g) J)® (cr^ «) tJ^) j . (17) 

One can show [t30|j5l| ] that frequency operators of the form (^) naturally imply quantum mechanical laws of large 
numbers (both weak and strong). 

Let us note that the sub-beam of N particles #1 (i.e. those going to the "left") is described by the product of 
reduced density matrices 

P{tj„tl}W ^ Ptltl{t)(^ ---(^ Pt^,t^{t) (18) 

which does not depend on times of detection of particles #2. This follows trivially from the fact that reduced 
density matrices of particles #1 do not depend on the form of Hamiltonian of particles #2, the fact whether the 
Hamiltonians are time-dependent or not being irrelevant. The switching-off functions are an element of definition of 
the Hamiltonians. 

Following exactly the same reasoning we can describe ensembles of n-tuples of particles and the associated n-time 
measurements. 



V. TWO-TIME MEASUREMENTS FOR ENTANGLED STATES IN NONLINEAR QUANTUM 

MECHANICS 



The above discussion does not really bring anything technically new if one sticks to linear evolutions of states. 

However, if nonlinear evolutions are involved the first two steps of our analysis, (^) and (^, are doubtful. The reason 
is that it is not clear whether the masurement performed on particle #1 should be just a simple projection producing 
a new initial condition at t = ti for particle #2. In linear quantum mechanics the projected state is an eigenvector. In 
nonlinear theory there are many inequivalent definitions of eigenvectors Moreover, if Doebner-Goldin equations 
are concerned the projection postulate should be formulated in a gauge independent way |^] 

Another point is that nonlinear dynamics is very sensitive to modifications of initial conditions and in the generic 
case the reduced density matrix of particle #2 will depend on the choice of both ti and the projector which performs 
the reduction. The well known Gedankenexperiment of Gisin is precisely a two-time measurement involving steps 
®"® ^^'^ ^ nonlinear evolution for ti < t < t2. It is known to predict an unphysical behavior of the particles and 
some kind of "faster-than- light" effect. 

One should know that an appropriately (or, rather, inappropriately) defined nonlinear dynamics of entangled states 
can produce similar unphysical phenomena jsj even if one does not explicitly use the projection postulate at the level 
of solutions. The Weinberg formalism j|] incorporates the projection postulate already at the level of Hamiltonian 
function of the composite system. In this way one defines a multiparticle nonlinear Schrodinger equation whose 
solutions imply probabilities one would have derived from the two step procedure (^-(^. 

It was shown later by Polchinski j^] and others |^,|lO| how to redefine the dynamics in a way which eliminates the 
unwanted properties of the Weinberg formalism. To get rid of the unphysical effect associated with the projection 
postulate Polchinski decided to eliminate the projection postulate itself. But then it turned out that strange effects 
occur at the level of the Many Worlds Interpretation and the "EPR phone" is replaced by an "Everett phone" . One 
cannot completely disagree with those who are not yet fully satisfied by the Polchinski formulation P,p| j4l| , ^ 

In the next section we propose a new framework, a sort of "golden middle" between Gisin and Polchinski. Although 
we keep the projection postulate and do not need Many Worlds (hence no Everett phone) we nevertheless maintain 
the local propeties of the Polchinski formalism (no EPR phone). 



VI. SWITCHING-OFF-FUNCTION MODIFICATION OF POLCHINSKI'S MULTIPARTICLE 

EXTENSION 

Nonlinear dynamics in the Wcinbcrg-Polchinski quantum mechanics is given by Schrodinger equations associated 
with non-bilinear Hamiltonian functions, 

Olpk 

It is essential that not all the possible Hamiltonian functions are acceptable. One accepts only those 7i(V', V') which 
can be written as 
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W(V'>)=7^(p)|^^l^,^^^,|. (20) 

For example 



is acceptable, whereas 



H(7A,Vi) = (V'+V'-+V5-^+)' (22) 

is not (as opposed to the acceptable functions the latter is not invariant under i-^ e*"|V'))- In linear quantum 
mechanics all Hamiltonian functions can be written as 

nii;,i:) = {mm - Trpi?|^^|^^^^, (23) 

and, hence, are acceptable. 

If we have two such particles each described by its own Hamiltonian function 

^i(V'i>^i)=^i(p)|,H^,)(^,(, (24) 

H2(V2,Vi2)=W2(p)|^^|^^>(^^|, (25) 

then the two-particle Hamiltonian function is simply their sum evaluated in appropriate one-particle states of particles 
#1 and #2, respectively. 

What makes this formulation ingeniously simple is the fact that for a generic entangled state 

l*) = E*feife^l^i)l^2) (26) 

fclfc2 

representing the two-particle system, states of the one particle subsystems may be represented by reduced density 
matrices 

Pi= J2 '^k,k.,^hkjk,){h\, (27) 

kilik2 

P2^ *'oife2*fel/Jfe)(^2|. (28) 

Due to the above mentioned acceptability condition it makes sense to consider 

Hi+2{^,^^HM\,,+n2{p)\^^. (29) 
The two-particle Schrodinger equation has again the Hamilton form 

^^^1^2 = — • 30) 

Having found its solution we can write with its help reduced density matrices of the subsystems. It can be shown in 
several different ways and at different levels of generality [p]J^ jri| , p^ that the dynamics of a reduced density matrix 
of one of those subsystems is independent of the choice of Hamiltonian function of the other subsystem. Therefore it 
is not posible to influence particle #1 by measurements performed on particle #2, and vice versa. This establishes 
locality of the formalism. 

Applying the above technique to the EPR situation we have no difficulty with describing measurements of, say, 
CTz ® (^z if the measurements are performed at the same time t = ti = t2. The formula for probabilities is the standard 
one: 

(*(t)|£;±®s±,|*(i)). (31) 

But what about two-time measurements? We do not have to our disposal the projection at ti since we know that it 
will produce unphysical nonlocal effects. 
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What we can do, however, is to resort to the switching-ofF function formulation which, as we have shown above, is 
equivalent to the one based on the projection postulate in the linear case. We define the two particle Hamiltonian 
function by 

Ht.Mit, *) = (^{t ~ + Sit - t2)H2{p)\^^. (32) 

The Schrodinger equation for two particles is again 

**fcife2 = • 33) 

It is important that, similarly to the linear case, the reduced density matrices of the subsystems depend only on their 
"own" tiS. This is a straightforward consequence of locality of the Polchinski formulation. 



VII. EXAMPLE: EVOLUTION OF A PAIR OF SPIN-1/2 PARTICLES 

The Hamiltonian functions in this example are 

Hi(V'i,V5i) = ^(V'ik.|V'i)V2 (34) 
n2{^2,^2) = B{il^2\a,\^2f/2. (35) 

The Polchinski two-particle extension is 

Hl+2(«',^') = ^(*|cr^®/|«')V2 + 5(1'|/®cr^|*)V2 (36) 
and the two-particle Schrodinger equation derived from this Hamiltonian function is 

zl^-) = (g)I + (g) a^l*)/ a^) |*). (37) 

Our modification is 

Ht,,t,(*, 'I') = eit - h)A(^\a, ® /|*) 2/2 + e{t - t2)B(^\I ® c7,|*) 2/2 (38) 

and 

i\if) = (^e{t - ti)A{^>\a, (g) I\'f)a, (g) I + 9{t - t2)B{'f\I (g) cr,\'f)I(g) a.'^l'i/), (39) 
with the general solution 

i (t)) ^ g-^^(*(0)l'^-®-f|*(0))'^-'»-f X5*e(T-tl)dT-lB(>I'(0)|7®<T,|*(0))7(»<T,/Je(T-t2)dT|^^Q-j^ 

— g-iA(CT^(0))l(T^K(t,ti) ^ g-i-B(o-^(0))2Cr^K(t,t2)|v[/(0)) (40) 

where {<Tz{0))k — Tr (^^(O)^^), K{t, tk) = /J 6{t — tk)dT. ( ^ ) describes the entire history of the two particles: From 
their "birth" at t = to their "deaths" at i = and t = t2- The standard Polchinski solution is recovered in the 
limits ti,t2 +00. 

There is absolutely no ambiguity in the switching-off-function formulation. An experimentalist who wants to know 
what are the theoretical predictions for his experiment has simply to look into his data and insert the detection times, 
ti and t2, into (^). 

An experiment with the beam consisting of N pairs is described by (^6|) with each of the entries of the form ( |40| ) 
but with initial conditions taken at t^s. 

To simplify further analysis we can assume that for all the pairs the times of flight At^ = t\ — tg, k — 1,2, 
i — 1, . . . ,iV, are the same and equal Atk- Under such assumptions averages of observables, say, X g)Y, can be 
computed as 

{Xg)Y) {'fAUAt.MX g> F|*Ati,At2(i))- (41) 

Averages of one-system observables, say X, are computed in the standard way 
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{X(E)I) 

= Tr (^e*'^<'^^("^^i'^'-''(*^^*i)Xe~'^^'^^(°)^^'^'-''(*''^*iVi(0))- (42) 

The average does not depend on At2. As we have aheady said this is a consequence of general local properties of the 
Polchinski extension. 

It is instructive to compare our result with the one we would have obtained on the basis of standard projection- 
postulate computation. Assume that until the first measurement the evolution is the same as before. At i = one 
performs a measurement of spin in some direction a, i.e. the observable is A" = a • cr. The corresponding projectors 
are E± — {I ± X)/2. Again we divide the calculation into steps. 

• At t = ti the state is 

\-^{ti)) = e-'^<'^-(o))i'^-*i Oe-''-^<'"^(°)>^''^*i|«'(0)). (43) 



Using the Zeno-type argument we reduce at t = ti by 



• We allow the "right" particle to evolve for ti < t < t2 but starting at ti with the initial condition ( |4^ ) 

\^±it2)) =/(g,e-^-«<*±(*i)l-f®'^-l*±(*i»^-(*^-*i)|^'±(ti)). (45) 

• Compute the conditional probability 

{^{tl)\E± (g) e*-^^*±(*l)l^®'^-l*±(*l)>'^-(*2-tl)£'_j_,g-iB(*±(ti)|J0o-,|'I'±(tl))cr,(t2-tl)|\J,(-^^)^) 



{^{ti)\E±<»I\^{ti)) 



(46) 



• The numerator of ( [46| ) is the joint probability we are looking for: 
where 

J (^2)) = (.-iM'i'm<^.(Si\-i'm'^.ti ^ g-iS(*±(ti)|/®<T,|*±(ti))<T,(t2-ti)g-iS(*(o)|/®<T,|*(o))fT,ti 1^(0)) (48) 

The latter formula can be compared with the one we have derived on the basis of our formalism and which can be 
written as 

|*ti (2(^2)) = e"*^^*^"^''^^'^-^'*^'')^'^^*! (8) e-^-S<*('')|7(»(T,|i'(o))<T,(t2-ti)g-i-B(i'(o)|/®CT,|i'(o))<T,ti|^(^Q^^^ ^^^g^ 

Figs. 2 and 3 show averages of ® ax (solid), ax® I (dashed), and I ® ax (dotted) calculated by means of the two 
procedures. The initial state is 

l*(0)) = ^|l)|2)-^|2)|l) (50) 



IDHlll'lV I2>=f-!i^/?). (51) 



where 

" sin(V8) ; ' ^^1- \ cos(V8) 

The parameters in Hamiltonians are A — 'H^ B = 1/2, and the detection times are ti = 3.5 and t2 — ^ (all in 
dimensionless units). In Fig. 2 we have used the switching-function approach. The dotted line representing the 
average oi I ® ax does not "notice" the measurement performed on particle #1. In Fig. 3 one can observe a slight 
change in the doted curve at t = ti. This is the nonlocal effect of the type described by Gisin. Until t = ti the 
evolution is described in the Polchinski way. As one can see the Zeno-type reasoning leads even in this case to the 
nonlocal influence between the two particles. The switching-function-modified Polchinski description is free from the 
difficulty. 
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VIII. PRE-SELECTION, POST-SELECTION, AND TELEPORTATION 



It seems we should finally make it clear how does our procedure affect the issue of teleportation of quantum states. 

Tclcportation is precisely an experiment of the type we have disciisscd: There arc correlated systems and information 
on results of measurements on one subsystem is used to create at later times an appropriate state of the other one. 
We claim there is no actual difficulty but one has to be careful with identification of pre- and post-selected ensembles 
if nonlinear evolution is involved. 

The simplest example of quantum teleportation is the following procedure. 

• A source produces a pair of particles in a singlet state. 

• An observer #1 ("Alice") performs a measurement of spin of particle #1. 

• If the result is +1/2 she phones to "Bob" and instructs him to remove particle #2. 

• If the result is —1/2 she phones to "Bob", and tells him to use his particle for further measurements. 

• As a result Bob obtains a beam of particles with spin +1/2. 

In more sophisticated versions of tclcportation one can use several particles and measurements on both sides may be 
accompanied by various additional operations. What is common to all the teleportation schemes is the presence of 
the "classical channel" (a phone, say) which allows Alice and Bob to exchange information and instructions. 

The above procedure should be contrasted with the correlation experiment in which Alice and Bob simply perform 
measurements and afterwards , when the experiment is completed, reject all pairs of data where spin +1/2 was found 
for particles #1. After having rejected a part of the data what is left may be regarded as an experiment performed 
by Bob on the sub-beam of particles #2 whose initial spin was +1/2. 

The teleportation scheme involves the so called pre-selection (Bob deals with a partial, or pre-selected, ensemble). 
Correlation experiment is based on post-selection (Bob first deals with the entire ensemble and post factum chooses 
the interesting part of the data). 

Linearity of quantum mechanics allows one to formally identify the two procedures in spite of the fact that they 
involve completely different experimental arrangements. To sec the differences one encounters in nonlinear situation 
consider the nonlinear Schrodinger equation of the type we have used in the example. 

Let the nonlinear Hamiltonian of particle #1 from our pair be 

A(*|cr (8) Jl^-) • cr (8) / ~ B(pi) • cr (8) J (52) 

where pi is the reduced density matrix of particle #1 and we have introduced the magnetic field 

B{pi) - TrpiCT. (53) 

which is proportional to the average magnetic moment of the beam. The Hamiltonian represents a mean-field inter- 
action of a single particle with the magnetic field created by the entire beam. 

In the singlet case the average magnetic field of the entire ensemble of particles #1 is zero. Alice may or may not 
perform her measurements but this is irreleveant unless she instructs Bob to undertake concrete filtering operations. 
However, if Bob pro-selects a half of the ensemble then obviously the ensemble will produce a nonvanishing magnetic 
field. The data collected with pre- and post-selected ensembles will be different. Each time Bob filters out a particle 
on the basis of Alice's information, he creates a new initial condition. Preparation of initial conditions must be 
accompanied by a projection. 

IX. DISCUSSION 

The procedure of computing multiple-time correlation functions we have introduced is not intended as a solution of 
the "measurement problem" . Moreover, there are other limitations of the procedure, especially arising in relativistic 
contexts. To give an example, it is not clear how to parametrize the switching-ofF functions for relativistic equations. 
However, from the very beginning our approach was not meant as a fundamental one. We have concentrated only on 
the part of the dynamics allowing the energy to leak out to detectors. In this respect we are doing neither worse not 
better than all the other approaches, both quantum and classical, that allow time dependent Hamiltonians. What we 
have tried to do was to clarify an important thinking error which seems to have plagued the literature on nonlinear 
Schrodinger equations. 
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There are many dangers one encounters while trying to extend Unear quantum mechanics to nonHnear domains. 
The reason is that Unear theories are, in a sense, very pathological. One can imagine the difficulties we would have 
in classical mechanics if the only potentials experimentally observed in Nature were those of harmonic oscillators. No 
doubt there would be some "impossibility theorems" about nonexistence of other potentials. 

This is precisely what happens in linear quantum mechanics since linear Schrodinger equation is mathematically 
equivalent to a classical Hamiltonian system describing uncountably many linear harmonic oscillators. Nonlinear 
Schrodinger equations correspond to non-quadratic Hamiltonian functions. 

Still, there is no way of stopping people use nonlinear Schrodinger equations. They are simply too useful and 
too interesting, and often very natural from a mathematical point of view. However, once one realizes that it may 
be justifed to contemplate nonlinearly evolving states, one immediately encounters another logical trap: A generic 
quantum state is not represented by a ray or a vector from a Hilbert space, but by a reduced density matrix which is 
not a one-dimensional projector. States represented by rays are a rarity enjoyed by completely isolated systems 
This is one of the most fundamental consequences of the existence of entangled states which, as we all know, are what 
makes quantum mechanics so special. 

It seems therefore that the psychological difficulty we have to face is that once we accept a possibility of nonlinearly 
evolving states, we have to immediately forget about any fundamental importance of nonlinear Schrodinger equations. 
This was realized by Polchinski and this is his great contribution to the subject. One has to start with nonlinear 
generators of evolution defined on the set of generic states, that is, density operators. But then the fundamental level 
of description must be the von Neumann one. 

There is nothing wrong with such a viewpoint. The set of solutions of nonlinear von Neumann equations is a Poisson 
manifold and, hence, a classical phase space with a Lie-Poisson dynamics. In this respect nonlinear von Neumann 
equations share many similarities with rigid-body or hydrodynamical equations, and even with the Nahm equations 
from Yang-Mills theories . Solutions of von Neumann equations are pure states (in the classical dynamical sense) 
on a Poisson manifold in the same sense as solutions of nonlinear Schrodinger equations are pure states on a Kahler 
manifold P,p"8|-p0t . Similarly to nonlinear Schrodinger equations the von Neumann ones (at least a surprisingly rich 
class of them) are integrable by Lax-pair and Darboux methods [^,|5| . 

It is clear that a density operator which solves a nonlinear von Neumann equation has the same ontological status 
as a state- vector solution of a Schrodinger equation. In our approach the density operator represents a state of a 
single quantum system and not of the entire ensemble. The ensemble is represented by a (finite or infinite) tensor 
product of single-particle states. In this respect we are not really very original since several other authors gave in 
different contexts (weak measurements [Q, quantum logic p7||) arguments for the same status of reduced density 
matrices (derived via reduction from entangled states) and state vectors. 

One dimensional projectors may be regarded as forming a kind of boundary of the Poisson manifold of pure states. 
To define a flow on this manifold it is not enough to define Hamiltonian functions on the boundary. Such a possibility 
indeed exists in linear quantum mechanics, but this is an example of pathological properties of linear Hamiltonian 
systems. But of course, once one has the dynamics on the entire manifold, one can contemplate its restriction to 
the boundary. And then one arrives again at nonlinear Schrodinger equations although somewhat stripped of a 
fundamental importance. This is actually what we have done in the examples discussed in this paper. 
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X. APPENDIX: NONLINEAR AVERAGES 

We define and use averages in the standard way: Averaging is linear and random variables are represented by 
Hermitian operators. In this way we avoid the arguments given by Mielnik p2| . But Hamiltonian functions are not 
represented by linear averages, otherwise the evolution would be linear. So to seriously face Mielnik's objections we 
have to say something about probability interpretation of nonlinear observables. 
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We are aware of at least two situations from a purely classical domain where several different ways of averaging 
(linear and nonlinear) may coexist simultaneously. What seems important, similarly to Hamiltonian equations of 
motion, they are related to optimization (i.e. variational) problems. The first class of examples is due to Renyi and 
is associated with his a-entropies. 



A. Kolmogorov-Nagumo averages and Renyi entropies 

Let E be the disjoint union of the sets Ei, . . . ,En having Ni, . . . , Nn elements respectively (X]fc=i ~ Let 
us suppose that we are interested only in knowing the subset Ek- The information characterizing an element of E 
consists of two parts: The first specifies the subset Ek containing this particular element and the second locates it 
within Ek- The amount of the second piece of information is, by Hartley formula, log^j iVj. On the other hand, 
to specify an element of E we need log^ N units of information. The amount necessary for the specification of the set 
Ek is therefore 

N 1 

h = log, N - log, Nk = log, — = log, — . (54) 

Nk Pk 

It follows that the amount of information received by learning that a single event of probability p took place equals 

I{p) = log, (55) 

This random variable is very interesting and unusual. It can be regarded as a classical example of a nonlinear 
observable. 

In statistical situations measured quantities correspond to averages of random variables. Therefore the average 
information is 

I = Y.Pk log,-. (56) 

k 

This is the Shannon formula and / is called the entropy of the probability distribution {pi, . . . ,p„} | |53[ |. 

Renyi gave examples of information theoretic questions where measures of information are those obtained by more 
general ways of averaging — the Kolmogorov-Nagumo (KN) function approach |Q . Let be a monotonic function. 
The KN average information can be defined by means of ip as 

If the generalized information measure is to satisfy the postulate of additivity, the function ip has to be linear or 
exponential (up to an affine transformation cp i— > aip + b which leaves all KN averages invariant). The linear function 
corresponds to Shannon's information. Choosing ip{x) — find 



^-'[i:Pk^{^og,^)y^lo^Y.P^). (58) 



Formula ( pq ) describes Renyi's a-entropy. 

The nonlinear Hamiltonian functions we have used in the example are of the form 

H^{p+-p-f 

-(p-i((^(9)p+ + (p(l)p_) -4(p+-p_)-4, (59) 

with p± = IV'iP: — V^- It follows that the Hamiltonian functions consist of two averages, one linear and one of 
the KN type. 

For general density matrices KN averages are 

{X(g)Y)KN = 'p'\Trp(p{X(g)Y)) (60) 

{X^I)kn^V'\T^i-PV{X<»I)) (61) 

^ (f-^ Ti- p(p{X) (g) I) (62) 

= ^-i(Tripi^(X)) (63) 

the latter being consistent with the Polchinski definition ( Tr i is the partial trace over the subsystem) . 
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B. g-averages and Tsallis entropies 



Tsallis non-extensive thermodynamics pS,B9] is based on the entropy 



and internal energy 



l-q 



H, = ^fi^^:Y,Pi'^^E,. (65) 



Thermodynamic equihbria are calculated by the standard variational principle which uses g-free energy Fq = Hq — TSq 

Jk'i 

normalized to unity. 

The definition of entropy becomes more natural if one recalls the formula 



and the constraint J2kPk — 1- Therefore there are two sets of probabilities in such a theory, and P^'^\ both 



p«-p = j3e(«-i)'"P-j3W {q-l)p\np. (66) 

This also sheds some light on the meaning of the parameter q in Tsallis' definition: q measures deviations from 
exponential and logarithmic functions. Thinking of the Shannon derivation of entropy we have given in the previous 
subsection one can see that q has something to do with possibilities of dividing ensembles into subensembles. That 
this is the case can be seen also in the formula for q-entropy of two independent systems {A and B) 

Sq{A UB)= Sq{A) + Sq{B) + (1 - q) S qiA) S q{B) . (67) 

Tsallis statistics has a wide range of applications. There are physical systems that are well described by g's quite 
far from unity. In plasma one finds q ~ 0.5, in multiparticle production process (hadronization) the best fit of 
experimental data yields q « 0.72. q-statistics arises whenever in the system one encounters long-range correlations, 
memory effects or fractal structures. 

In a sense even more interesting from our perspective are purely classical applications of q-statistics since they 
provide many insights into the meaning of non-linear averages. Among such applications one finds statistics of goals 
in football championships, citations o scientific papers, the Zipf-Mandelbrot law in linguistics and sociology, and many 
others. 

From the point of view of nonlinear quantum dynamics it is interesting to look at neighborhoods of Tsallis ther- 
modynamical equilibria. A variational principle based on q-averaged energy leads to the nonlinear von Neumann 



equation |37 



ip^[H,p'^] (68) 
and free energy Fq plays a role of a stability function which characterizes a dynamical equilibrium. 
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Detectors 

FIG. 1. A history. Projector Ei does not "reduce" states. 
Filters Ui and U2 are represented by unitary operators. 
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FIG. 2. Averages of the three obscrvables in the switch- 
ing-function formulation. The dotted line shows the evolution 
of observable ax associated with particle #2 which is detected 
at t = t2 = 8. Earlier detection of particle #1 at ti = 3.5 
does not influence particle #2. 
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3. Averages of the three observables calculated in 



Zeno-type way. Measurement at t = ti = 3.5 performed on 
particle #1 nonlocally influences the behavior of particle #2. 
As opposed to the plot from Fig. 2 the dotted line is modifled 
at t = 3.5. This is Gisin-type nonlocality. 
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